res_table <- xtable(res3[,c("label", "coef", "p.val","p2", "treat")])
names(res_table) <- c("Variable", "Estimate", "Uncorrected p-value", "Benjamini-Hochberg corrected p-value")
print(res_table, include.rownames=FALSE, file= "results/tablee2.tex")
